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Oh. The distribution of W and Z bosons produced with small transverse momentum (p^) 

^ ' at hadron colliders receives important contributions from large logarithms arising from 

soft gluon emission. Although conventionally the all-orders resummation of these 'Su- 
dakov' logarithms is performed in impact parameter (Fourier transform) space, pT-space 
. resummation is also possible, and offers certain advantages. We present a detailed phe- 

$H ' nomenological analysis of W and Z production at small Pt at the Tevatron pp collider, 

^ using pT-space resummation. A good description of the CDF and DO data can be ob- 

tained provided a significant non-perturbative contribution is included. We also present 
predictions for the LHC. 
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1 Introduction 



The description of gauge boson production at hadron colliders has recently attracted much 
theoretical interest, especially in the light of future high precision experiments at the Tevatron 
and the LHC |l|, ^, . Reliable predictions can only be obtained if soft gluon radiation effects 
are correctly taken into account. Theoretically the soft gluon emission manifests itself in the 
presence of large logarithmic corrections (Sudakov logarithms). For the particular case of 
the transverse momentum (p^) distribution of a boson produced with invariant mass Q, 
the Sudakov logarithms are the logarithms of the ratio Q"^ /p%- In the small p-^ limit the 
logarithms diverge and the standard fixed-order perturbation theory approach breaks down. 
However, a finite result can be recovered if the soft gluon emission is accounted for to all 
orders in as- This is achieved by resumming the logarithmic corrections. 

Resummation can be performed either directly in transverse momentum (pt) space or 
in the Fourier conjugate impact parameter (6) space. The most leading logarithmic contri- 
butions, of the form a^\r?"'~^{Q'^ /p\), can be directly resummed in px space (the so-called 
Double Leading Logarithm Approximation) Q. The impact parameter (6) space method 
allows one to resum subleading logarithms, including those 'kinematic' logarithms arising as 
a direct result of transverse momentum conservation @. Although very successful theoret- 
ically, the b space method suffers from certain deficiencies and drawbacks which need to be 
'fixed' in order to obtain a satisfactory agreement between the theoretical predictions and 
experimental data. For example, one experiences difficulties when matching the resummed 
(small Pt) and fixed-order (large px) predictions. Moreover, it is impossible to make predic- 
tions for any value of without a prescription dealing with the non-perturbative regime of 
large These difficulties can be naturally circumvented if the resummation is performed in 
Px space Unfortunately the px space methods which have been developed so far for re- 
summing subleading logarithms are derived from the b space approach |9|, |l^ , and as such 
they simply provide an approximation to the b space result. However, the goal is to develop 
a phenomenologically useful px space expression that reproduces all the good features of the 
b space resummation without the drawbacks related to this method. 

In a previous paper ||l^ we proposed a px space resummation formalism at the parton 
level which resums the first four 'towers' of logarithms (i.e. terms of the form a" \t?^~'^ [Q"^ / p\) 
m = 1, ..,4), including the effects of transverse momentum conservation. The differences be- 
tween our formalism (KS) and other px space approaches (FNR [^] and EV [Q) were discussed 
in [11|. Here we want to concentrate on the practical applications of our formalism, in par- 
ticular on the comparison with available Tevatron data. 

In our analysis we use the most recent sets of CDF [12| and DO data [13| on Z production 
and DO data |14| on production. Since we are only interested here in the resummed part of 
the cross section and do not perform matching with the fixed-order part, we do not consider 
data above px > 25 GeV. In this px range the resummed part accounts for almost the entire 
cross section. 

In a manner similar to the b space formalism, the px space formalism is incomplete with- 
out a prescription for dealing with the non-perturbative effects. Indeed previous phenomeno- 



^ Recently it has been argued that the b space drawbacks discussed here are less relevant at Tevatron 
and LHC energies. 
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logical analyses have shown that the very small Pt region is dominated by non-perturbative 
contributions. Here we use the method of introducing non-perturbative effects in space 
first proposed in [^. We investigate the form and size of the non-perturbative contributions 
obtained from fits to the data. We finish our investigations by commenting on W and Z 
boson production at the LHC. 



2 Theoretical cross section for pp ^ Z + X 



In this section we summarise the derivation of the main theoretical results for space used 
in fits to the data. For the discussion of the b space results the reader is referred to [l5| . 

The resummed part of the theoretical cross section in space for a Drell-Yan-type 
process follows from the b space formula, cf. [0] 



da 



^ q 



1 

2 Jo 



dxA dxB S { XAXB 



dbbJoip^b) exp[5(6,Q)] /g/^ (^^'t) [^^ 



bo 



bo 



with o"o = 47ra^/(9s), bo = 2exp(— 7^), and where 



S{b, Q2 



dfi^ 



\n[%]A{as{pi')) + B{as{-^')) 



2tt 



as 



27r 



(2) 
(3) 



i=l ^ ' 1=1 

We first consider the non-singlet (NS) cross-section, i.e. we introduce 

/(j/H — fq/H ~ fq/H 

as modified higher-order NS parton distributions. The modified parton distributions are 
related to the MS parton distributions, /, by a convolution 16, 15 1 

fa/ nix A, ^J) =^ I ^Cac(^,fj) fc/H{z,f^) , (4) 

where (a, b ^ g) 

Cag{z,fJ-) = as{fJ,)Tji 2z{l - z) 



1-^+ y-4U(l-z) 



and asifi) = ^,Cf = 4/3, Tr = 1/2. 
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The A^-th moment of the cross section with respect to r = Q^/s has the form 



M{N) 



dr T 



el\ dbbUp^b) exp [5(6, Q)] /'^/^{N, ^) 4s(iV, ^) 



(5) 



Solving the DGLAP equation for the A^-th moment of the modified parton distribution 
f'g/ni^^Q) ~ J dxHX^f'qiu{xH,Q), and integration by parts lead to (cf. 



M{N) 



dp^ 



Jo 



dx Ji{x) exp 



Six,Q)-2 



7jv(as(M)) 



where x = p^b. 

In order to obtain an expression for the hadron level cross section, the following approx- 
imation is introduced 



exp 



[Pt dff 

Six,Q) -2 y,2p2 -32-7^(as(M)) 



exp [S{x, Q)] . 



(6) 



The above equality is exact for the first four towers of logarithms; it is only the fifth 
tower that contains the first modified anomalous dimension coefficient 7^^'*- This can be 
easily seen by expanding the exponential in @ (assuming here a fixed coupling constant for 
simplicity) 



exp 



E 

Af=0 



(-1) 



N 



liA^^^ds + A^^^ds^ + ...){L + Lbf + (S«a-, + B^^^ds^ + ...)(L + Lb) 



+2{j'f,'^ds+j';,^'ds-' + ...)L, 



'(2).- 2 



N 



where L = ln(Q^/p^) and = ln(x^/6Q). The first term containing 7^^^ which does not 

vanish after integration over x is of the form q^^^''^^ same statement 

holds also for the singlet parton distribution functions. 
The resulting expression 

d 



M{N) = ^1 f'^,j,{N,p^)~f'g/j,{N,p^) dxJi{x) exp[5(x,g)] 



(7) 
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can now be transformed back to momentum space by the means of the inverse MeUin trans- 
form 

da (To ^ 2 1 re/ 



dp], dQ"^ 



Jt\ dxJiix)exp[S{x,Q)]fg/^{xA,PT)f'q/BixB,PT)\- (8) 



dplj. 

At the parton level we calculated the quantity^ 

1 d 

= -^exp(5^(Q)) dxxJi{x)—exp[S{x,Q)] 



"T" N=l V " / m=0 "- k=0 

N-m-k-1 -I N-m-k-l-l N-m-k-l-j~l 

X 



-m— fe— 1 -j^ N~m—k—l—l ^ N —m—k—l 

X! 77 77 X! 



, ^ /! ^ j! ^ i!(iV-m-A:-/-j-z-l)! 
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^ m k I ^ i N— m—k—l— 

X C2 C^C^Ct^CqCi nCnT]\f^.in+2k+3l+Aj+5i+n-2 

n=l 

= (9) 
in terms of resummed towers of logarithms in space. Here S{x,Q) = S{x,Q) — Sr]{Q) 

roci 

with Srj and c coefficients defined in The expression for / dxJi{x)exp[S{x,Q)] can 
be derived in a similar manner 

/ dxJi(x)exp[5(x,Q)] = exp(5r,) 2^ 2^ — T 2^ 77 

N-m~k-l N-m-k-l-1 N-rn-k-l-j-l ^ 



w m k I J i N— m—k—l— j—i—1 

X C2 C^C^C^CqC^ TN+m+2k+3l+Aj+5i-l 

= s2(pT,g), (10) 

Finally we arrive at the p^ space formula for the Drell-Yan cross section at the hadron level 
da 



dp], dQ'^ 



{^2{Vt. Q) f'q/A{^ A, Pt) fq/Bi^B^PT)] ■ (11) 



^see Eq. (21) in |1C 
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In what follows we will refer to the result in Eq. (|Il|) as the KS hadron-level formula in 
space. 

In principle, the parton level formula (^) allows us to resum any number of towers of 
logarithms. In practice, however, the fifth tower of logarithms cannot be fully taken into 
account due to the lack of knowledge of the coefficient A'^'^\ Since our approximation (|6|) is 
valid only up to the fifth tower too, Eq. (pT]) can be used to resum the first four towers of 



logarithms, in other words the summation in (|lO|) stops at = 4. In |10|, the contributions 
from fifth and higher towers were estimated to be numerically very small in the region of 
of interest. 

The analogous expression for the transverse momentum distribution of a massive vector 
boson V produced va. pp + X \s 



'T 



qq' 



d 



{S2(Pt, My) /^/^(xa,^5t)/;/s(xb,Pt)} , (12) 



where 



dp-r 



N 



rrV _ i M V = W^, 

where Vqq' denotes the appropriate CKM matrix element, and Vg, Ag are the vector and axial 
couplings of the Z boson to quarks. 

In practice it is convenient to split the differentiation in (|l^ ) into two terms 



T— = (^Oy^ U^n' [ dXAdXB^ \ XAXB - 

dpT 



qq' \ / 

- — ^l{pT,Mv)fg,A{xA,PT)fg>/BixB,PT) 
Pt 

+ T.2{Pt,Mv)-^ f'glA{XA,PT) f'g:iB{XB,PT)\ \ ■ (14) 



'T 



This trick allows us to apply an inevitable numerical derivative only to the product of the 
parton distributions and not to the whole expression, leading to a reduction of the numerical 
error. 



2.1 Inclusion of the non-perturbative effects in space 

The form of the non-perturbative ansatz in space is expected to be important only in 
regions where perturbation theory fails, i.e. at the very low values oi p-^ < 2 — 3 GeV. In 
contrast, the higher Pt region can be described purely by the resummed perturbative QCD 



5 



expression. We choose to incorporate the low energy effects using the form of the space 
non-perturbative function F^^{p'y) advocated in 



F''^{p,) = l-exp [-api]. 



(15) 



The role of this function is to account for the distribution in the very low px region, and here 
we are assuming that the shape is approximately gaussian. However in order to combine this 
with the perturbative result, the latter needs to be 'frozen' or 'switched off' at some critical 
value of where the coupling becomes large. A similar freezing is required in the b space 
approach where the coupling is effectively as{l/b). In other words we require not only (i) a 
form F^^{prj.) for the distribution in the non-perturbative region, but also (ii) a prescription 
for moving smoothly from the perturbative to the non-perturbative region. One possibility 
for the latter is the 'freezing' prescription of |S] 



IpI 



exp 



pI 



which has the property 



Pt* 




Pt lim 



Pt > PTlin 

Pt <■ Pthu 



(16) 



(17) 



It is important to note that there are two pieces of information contained in this definition: 
the value of the limiting value PTUm and the abruptness of the transition to this value. The 



use of a gaussian function in the definition (16), compared to say a power law function, 
implies a rapid transition that, as we shall see below, is consistent with the data. 
Applying the above prescription to our expression ( [T^ ) leads to 



T~~ = '^0^2 ^w' I dxAdxsS ( XAXB 
dpT ^ \ 



1 dpT 
Pt* dp^ 



-T.i{prr^,Mv)f'g/A{xA,PT*) fq'/BixB,PT*)F {p 



+ S2(Pt*, My) 



dpT* d 



dpT dpT* 

+ S2 (Pt* , My ) f^/^ (x A , Pt* ) /g//B (XB , Pt: 



fq/Ai^A,PT*) fq'/Bi^B,PT* 

d 



F'^'iPT) 



'■pNP 



dpT 



{Pt 



(18) 



Note that the simple form of F^^{pt^) in the present framework does not take into 
account a possible dependence on Q and x. This is in contrast to the b space treatment 



of [21 1, where an x-dependent linear term in b was added to the argument of the gaussian 
non-perturbative function. It has also been argued that the width of the non-perturbative 
gaussian distribution, in our case the parameter 1/a, should increase linearly with logQ 0. 

^ Indeed, fits of a non-perturbative gaussian distribution to the low energy data typically give much 
larger values of 5 (« 0.55 GeV~^), which suggests a strong dependence of a on Q. 
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Since in the present case we are only interested in W, Z production at a single collider energy, 
the values of Q and x are essentially fixed at My and My / \/s respectively. Therefore we are 
not able to say anything about the form of the dependence of the non-perturbative parameters 
on Q and x. Nor will we investigate different functional forms for F^^(pt) — the simple 
gaussian form in (|T^ ) will allow perfectly good fits to the Tevatron data. 

The lack of information on the x dependence of the non-perturbative contributions should 
be borne in mind when considering the predictions for W and Z production at the LHC. 



3 Results and discussion 

For the parton level cross section we have advocated [10| the use of the renormalization scale 



/x/j = p^^^'^ Q^^/'^^ as a means of eliminating certain logarithmic terms from the Sudakov factor 
and thus increasing the reliability of our approach. Since the renormalization scale determines 
the strength of the coupling in the theoretical predictions, it must somehow depend on the 
size of the transverse momentum and we require the choice of the scale to reflect this fact. 
Moreover, for values of px where perturbative QCD can be safely applied and for the values 
of Q considered here, such a is always bigger then the h quark mass, thus lessening the 
relevance of the correction due to the treatment of quark mass thresholds. Another obvious 
choice for the renormalization scale is = However, since we find only a very small 
dependence of the resummed part of the cross section on the choice of ///j, from now on we 
use = p^^'^'^ Q^^/^^ as the default choice for the KS approach. 

The Drell-Yan cross section (|Tl| ) has been derived in the limit of a fixed number of 
quark flavours, Nf, which implies that no quark mass threshold effects are considered. In 
the original h space approach, the dependence on Nj enters in the Sudakov factor through 
the coefficients and through the (3 function in the expansion of a^. For our 

space method we propose to change Nf according to the number of flavours active at the 
renormalization scale at which is calculated while the remainder of the expression is derived 
in the massless quark limit. With the choices of the scale we use this roughly corresponds 
to the energy scale of the emitted gluons and fits comfortably into the physical picture of 
the process. Changing the number of active quark flavours Nf as p-^ varies immediately 
leads to the problem of obtaining reliable predictions free of unphysical discontinuities. To 
overcome this we use an analytically extended scheme which incorporates finite-mass 



quark threshold effects into the running of the coupling, as proposed by Brodsky et al. |19] 



By connecting the coupling directly to the analytic and physically-defined effective charge 



scheme, the authors of |19] obtain an analytic expression for the effective number of flavours 



which is a continuous function of the renormalization scale and the quark masses. 



The results presented below are for the Tevatron experiments, CDF and DO, at -^i = 
1.8 TeV. Unless stated otherwise we use the factorization scale jif = p^, MRST99 par- 
ton distribution functions (central gluon), branching ratios BR{Z e^e^) = 3.366%, 

*Other choices of /x_r have been considered in the literature, for example the authors of jl^] proposed to 
take hr = \/p| + 



7 



BR{W — > ev) = 11.1% and the world average value of the strong coupling as{Mz) = 0.1175. 
To normalize the theory predictions to the data we take only those experimental points with 
< 15 GeV. 

There is a significant amount of Tevatron data on W and Z production that should, in 
principle, allow a precise determination of the non-perturbative parameters from fits to the 
data. However since the measurement of the W transverse momentum requires correcting 
for detector effects that are much stronger that in the Z measurement case, for the purpose 
of this analysis we take only the Z data. Again, we consider only those experimental points 
with < 15 GeV for the fit range. The overall normalization is taken as a free parameter, 
since we are primarily concerned with the shape of the distributions. |^ 

Generally we find that the best x^/d.o.f. value is obtained by values of \/l/a and pTUm of 
order 3 — 4 GeV. In this context the values proposed by the EV collaboration a = 0.1 GeV~^, 
Pxiim = 4 GeV provide one of the best fits and describe the Z data well. This is also in 
agreement with the CDF and DO analysis, cf. |12, 13]. Furthermore we find that there is a 



wide range of strongly correlated values of larger a and smaller pTUm for which x^/d.o.f. is 
only minimally worse. This is illustrated in Fig. |l|, which shows the equal contours in the 
plane for x^/d.o.f. = 1 and 0.75. 



- 0.6 

- 0.5 

- 0.4 



- 0.2 

I I I I I I I I I 

1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 



p T lim 



Figure 1: The contours of equal iii the a, pxiim plane for the KS space approach with 
the non-perturbative input of the form ( |T5| , |T^ ). Both CDF and DO data (with separate 
normalization) for < 15 GeV are used in the fit. The outer and inner contours correspond 
to x^/d.o.f. = 1 and 0.75 respectively. 

^The total W and Z cross sections are known to be well described by NNLO perturbation theory, see e.g. 

io). 
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The correlation between a and pTUm is easy to understand. Increasing pTUm corresponds 
to requiring the non-perturbative contribution to describe the data out to a larger value of px, 
and therefore a broader gaussian distribution (equivalently, smaller a) is required. The fact 
that the fit deteriorates sharply as pxiim is made very small shows that ('frozen') perturbation 
theory alone cannot describe the data over the whole range. 

Given the large variation in the allowed values of a and PTUm, it is difficult to gauge the 
predictive power of these results, especially when one allows for a possible additional x and 
Q dependence. We also find that with the current experimental data there is no need to 
introduce additional overall smearing, as proposed in [^]. A modification of -F^^(pt)) such 
as adding a linear term in the exponential or using a different freezing method, does not 
significantly improve the fit either. 

In Figs. ^, ^, ^ we present a comparison between experimental data on Z production 
as measured by CDF and DO, W production as measured by DO, and various theoretical 
distributions calculated using (i) the b space method, (ii) the EV p^ space method, and (iii) 
the KS Pt space method. We observe good agreement between the data and the theoretical 
predictions for all three methods, in the range of = ~ 25 GeV. In general, the b space 
distribution is more 'peaked' than the p^ space equivalents. This effect is, however, very 
susceptible to the choice of the non-perturbative function and values of the non-perturbative 
parameters. The b space distribution is also higher in the intermediate range of px = 10 ~ 20 
GeV, where the non-perturbative physics does not influence the resummed perturbative re- 
sult. In this region the KS distribution approximates the b space result better than the 
corresponding EV distribution. Given that the KS formalism resums more towers of loga- 
rithms than the EV formalism, this is an expected result. The increase of the cross section 
due to incorporating the fourth, NNNL, tower can be as big as 4% for some values of p-r, 



both for W and Z production |10|. Interestingly we also observe a significant sensitivity to 
the value of as{Mz) used in the calculations. A variation of as{Mz) by ±0.005 around its 
average value, 0.1175, can cause, for some values of p-j-, a more than ±8% change in the Z 
Pt distribution. 

The transverse momentum distribution of W^s and Z's at the LHC, predicted in the p^ 
space formalism (KS), is shown in Fig. ^. The results agree with similar analyses performed 
using the b space method @. For the sake of this analysis we used the standard Tevatron 
values of the non-perturbative parameters o = 0.1 GeV~^ and Ptu^ = 4 GeV in p^ space. 
This may prove to be a very unwise assumption if the non-perturbative parameterization 
does depend significantly on the partons momentum fractions, for example in the way it was 
proposed for the b space method ( [pl| ). However it does provide a useful benchmark and a 
reasonable 'first guess'. 

4 Summary 

We have applied the KS resummation technique inpx space, developed in ||l^, to the hadronic 
production of vector bosons. At the hadron level our approach retains the potential of the 
full resummation of the first four towers of logarithms. We also allow for a non-perturbative 
contribution, with a smooth interpolation between the perturbative and non-perturbative 
regimes at small p^- Our numerical results generally show good agreement with recent data 
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Figure 2: Comparison between CDF data on Z production and theoretical predictions for 
the b space method, space method in the EV approach and in the KS approach. For the 
b space method we use an effective gaussian form of as in |15]. 
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Figure 3: Comparison between DO data on Z production and theoretical predictions for the 
b space method, space method in the EV approach and in the KS approach. For the b 
space method we use an effective gaussian form of F^^ as in [15|. 
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Figure 4: Comparison between DO data on W production and theoretical predictions for the 
b space method, space method in the EV approach and in the KS approach. For the b 
space method we use an effective gaussian form of F^^ as in [15|. 
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Figure 5: px space predictions (KS formalism) for the transverse momentum distribution at 
the LHC of: (a) W boson, (b) Z boson. 

on W and Z boson production from the Tevatron coUider. 

For the resummed part of the da/dp-r distribution we observe rather weak dependence 
on the renormahzation scale, but some sensitivity to the value of as{Mz)- However the 
non-perturbative contribution, which at present must be determined from fits to data, is not 
well determined. We have not attempted to estimate the full error on the non-perturbative 
contribution. Rather, we showed that a simple Gaussian form, with reasonable values of the 
parameters, gives an acceptable fit. The width of the gaussian and the transition point be- 
tween the perturbative and non-perturbative regions are, however, strongly correlated. This, 
combined with the lack of knowledge of the x and Q dependence of the non-perturbative 
parameters, makes it difficult to formulate precise predictions for the corresponding distri- 
butions at the LHC. 
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